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SUMMARY 

A method utilizing source singularities is presented for obtaining 
the linearized downwash field due to lifting wings of infinitesimal 
thickness at subsonic and supersonic speeds. The distribution function 
for the source singularities is specified by the loading on the wing. 
The method is applied to derive generalized formulas for the downwash 
field due to uniformly loaded swept and rectangular wings at subsonic 
and supersonic speeds. The utilization of these formulas to obtain the 
downwash due to wings of arbitrary loading is indicated. An example of 
the procedure is given in which specific formulas are derived for the 
downwash field due to a rectangular wing at supersonic speeds for a 
uniform loading and for a linear chordwise variation in loading. 


INTRODUCTION 

Several methods based on linearized theory are available to obtain 
the downwash field due to lifting surfaces at subsonic and supersonic 
speeds for use, as an example, in stability calculations. The calcu- 
lation of the downwash field at subsonic speeds has relied almost 
exclusively on Prandtl's lifting-line theory, which is based on the 
concept of a horseshoe vortex (for example, reference 1) . Present 
methods for calculating the downwash field at supersonic speeds are 
those utilizing conical flows (reference 2) , potential doublets (refer- 
ence 3), vortices (references 4 and 5), and pressure doublets (refer- 
ences 6 and 7). The integrations required of the foregoing vortex 
or doublet singularities or in the coni cal -flow method in order to 
obtain exact solutions of the linearized equations for lifting surfaces 
have generally been found to be difficult} therefore practice has 
usually had recourse to approximate methods based on lifting-line 
theories (references 5 and 7). 

The present report, prepared at the NACA Lewis laboratory, indi- 
cates a method that is intended to facilitate the computations for 
obtaining the exact linearized downwash field due to lifting surfaces 
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at subsonic and supersonic speeds. The method utilizes source singu- 
larities with the distribution function specified by the loading on 
the surface. The method is applied herein to derive formulas for the 
dovnwash field due to uniformly loaded swept and rectangular wings of 
infinitesimal thickness at subsonic and supersonic speeds. The utili- 
zation of these formulas to obtain the downwash field due to wings of 
arbitrary loading, by means of the correspondence relations presented 
in reference 8, is then indicated. An example of the procedure is 
given in which specific formulas are derived for the rectangular wing 
at supersonic speeds for a uniform loading and for a linear chordwise 
variation in loading. 


SYMBOLS 

The following symbols are used in this report: 


A,B,C^,C 2 refer to regions bounded by foremost Mach aft cone 

p p ^ corners of rectangular wing at supersonic speed 

5 (plan view of regions in fig. 3) 


a 


B 

*L- 


constant used to describe prescribed linear chordwise 
variation in loading 

\ 

= ^M^ - 1 (also used to refer to Mach cone region as 

indicated in fig. 3) 

Nl - M 2 


b 


arbitrary constant 


c 


chord for rectangular wing 


h wing semispan 

refers to integral arid solution expressed by equa- 
tions (25b) and ( 25e) respectively 

I 2 refers to integral and solution expressed by equa- 

tions (25c) and (25g), respectively 

K constant in equation (1) (K * — at subsonic speeds; 

Git. 

K = — at supersonic speeds) 
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M 


free-stream Mach number 


m 


r 

r* 


cot A > where A is angle of sweep of wing leading 
edge (fig. 1) 



u,v,w 


disturbance velocities of fluid in flow field alpng 
x-, y-, and z-axes, respectively (fig. 1) 


u 0 


value of u on upper surface of wing for uniform 
prescribed loading 


refers to downwash due to left half of wing 


w, 


0 


refers to downwash due to uniform prescribed loading 
on wing 


refers to downwash due to linear chordwise variation 
in prescribed loading on wing 


w L contribution to downwash of continuous portion of. 

leading edge for semi-infinite oblique wing with 
uniform loading shown' in figure 1; also represents 
downwash due to semi-infinite line source originating 
at leading edge of center section 

v^j) contribution to downwash of discontinuity of leading 

edge at origin in semi-infinite oblique wing with 
uniform loading shown in figure 1 

w n ^ k downwash due to term in series formed by expressing u 

* on wing as function of x, y, and u 0 


contribution to downwash of streamwise side edge for 
semi-infinite oblique wing with uniform loading 
shown in figure lj also represents downwash due to 
semi-infinite streamwise line source originating 
at leading edge of center section 


w *o»yo 


represents downwash due to semi-infinite line source 
originating at (x Q ,y 0 ) on right and left halves of 

wing, respectively 
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rectangular coordinates with origin at leading edge 
of center section (fig. 1) 

oblique coordinates related to rectangular coordinates 
according to equations (5) 

rectangular coordinates indicating origin of semi- 
infinite line source 

y a = y~ h 

% \ 

y-b = y +h 

(Av ) contribution to dovnwash of wing cut-off at trailing 

® ^ edge for uniform loading 

e infinitesimal distance in y-direction across side 

edge 

f ' auxiliary variables used to replace x, y, z, x' , 

y ' , and z "respectively 

^ a - B 

^ b B 

£ d = V|B 2 |(y 2 +z 2 ) 

v d = VI B 2 |(y ,2+z ' 2 ) 

i 2 ,^2 upper limit for integral in equation (7) (See dis- 

c " cussion following equation (7).) 

T region of integration (See discussion following 

equations (l) and (3) .) . 

Q function representing solution of linearized partial 

differential equation (2) 

Subscripts : 

A,B,C, ,C 9 refer to corresponding regions indicated in figure 3 

j_ L, 

D 1' D 2> , ** D 5 



x,y,z 

x',y',z’ 

Wo 
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c r refers to chord at center section of swept wing 

c^ refers to chord at tip section of swept wing 

L,S,T,LD,TD refer to continuous portion of leading edge, side 

edge, trailing edge, discontinuity in leading edge, 
and discontinuity in trailing edge, respectively 

u,Z refer to upper and lower surfaces of airfoil, 

respectively 

Single or successive subscript coordinates indicate partial 
differentiation with respect to subscript variable. 

x s , yk indicates partial differentiation with respect to x 

and y, s and t times, respectively 


BASIC THEORY 

The analysis is based on the usual assumptions for thin airfoils 
in the linearized potential field. A solution for the disturbance 
parameters can thus be obtained by integrations of source and doublet 
singularities in the plane of the wing (z =0). The basic equation is 
(for example, references 9 and 10) 


^(x,y,z) 



( 1 ) 


where the function Q is a solution of the linearized partial differ- 
ential equation for subsonic and supersonic flows 



In equation (l), r = V(x4) 2 - B 2 [(y-T)) 2 +z 2 3 and the region T 
includes the entire z = 0 plane that can influence the point (x,y, z) . 
At subsonic speeds, the factor K is equal to l/2ir. At supersonic 
speeds, the factor K is equal to l/n and only the finite part of 
the integral is used. It is important to note that the function 0 can 
represent either the velocity potential or any of its derivatives, ar.d 
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if all these values vanish sufficiently fax ahead of the wing, the 
integrals of 0 are solutions of equation ( 2) . . In case the deriva- 
tive of ft becomes infinite at one or more points, the substitution 
of this derivative for 0 in equation (2) depends on the condition 
that the isolation of each singularity yields a finite integrand in 
the limit. 


SOURCE DISTRIBUTION FOR DOWNWASH FIELD 

A lifting wing of infinitesimal thickness is considered. The 
nerturbation velocity u vanishes everywhere in the z = 0 plane 
except on the wing itself, and the perturbation velocity v vanishes 
everywhere in the z = 0 plane except on the wing and in the wake. 

In equation (l) , letting ft = w and noting that 


u 


- w. 


= 0 


w z a B u - V|J 


<V« ‘ ' ( V* 


(T Ju - - <Vi 


result in 


w(x,y,z) 



(B 2 u^-v r( ) dfc dtj 


(3) 


where T represents one surface of the wing and of the wake that can 
influence the point x,y,z. The quantity v^ is a function of u 
through the irrotationality relation 



Therefore , 
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t i ■ S L u ' 1 ' 1 <4b) 

According to equation (3), the perturbation velocity w at any 
point in the flow field due to a lifting- surface is determined by an 
integration of elementary source solutions with the distribution 
function given by (B 2 U£-v^) . 

c 

THE SEMI -INFINITE OBLIQUE WING WITH UNIFORM LOADING 

An important application of equation (3) can be obtained by con- 
sidering a semi-infinite oblique wing with uniform loading (fig. l) . 

The origin of the coordinate system is taken at the intersection of 
the leading and side edges, where both edges extend to infinity. For 
uniform loading, the term u - in equation (3) vanishes everywhere over 
the region t except across the leading edge. In evaluating the disr 
continuities in u or v that occur across the edges, a limiting- 
procedure is used throughout the present analysis, which corresponds 
to the assumption of a linear variation in u or v across an 
infinitesimal strip of the edge. The distribution of v over the 
region t for this type of loading, as obtained from equation (4a) , 
is shown in figure 2, and v^ is seen to vanish everywhere except 
across the edges. In terms of equation (3), therefore, the downwash 
field for a semi-infinite oblique wing of uniform loading is obtained 
by means of a line integration of sources along the edges. 

It is subsequently shown that the downwash solutions for the semi- 
infinite line sources along the leading and side edges may be used by 
simple manipulation to obtain the downwash field for finite plan forms . 
of uniform loading. 


The Semi -Infinite Oblique Leading Edge 

The semi-infinite oblique leading edge with origin at (0,0,0) con- 
tributes to u| and v^ along the edge and, by virtue of the dis- 
continuity of the edge at the origin, contributes to v^ along the 
x-axis beginning at the origin and extending backwards to infinity. 

(See distribution of v in fig. 2.) The downwash contribution of 
the continuous portion of the leading edge is designated w^j whereas 
the contribution of the edge discontinuity at the origin is designated 

W LD* 
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Contribution of continuous portion of leading edge . - The integra- 
tion of equation (3) along the oblique leading edge is most conveniently 
performed in terms of an oblique system of coordinates, such that 


x' = x-mB 2 y 
y'= y-mx 

z’ - z A/L-m 2 B 2 


(5) 


It may be shown that the differential equation (2) is invariant under 
the change of variables x— *x' , y— ►y*, and z— >z' at both subsonic 
and supersonic speeds. If Q (x,y, z) is a solution of equation (2), 
^(x',y'}Z') is therefore also a solution. For other examples of the 
use of the oblique transformation in wing-theory problems, see 'refer- 
ences 11 to 13. 

In the oblique coordinate system, equation (3) evaluated along 
the wing leading edge becomes 


where 


Y^xS-y',*') « -K J* p-w z ,Ari , d^' 


r' = /^(x'-£') 2 - B 2 (y ,2 +z' 2 ) 


( 6 ) 


Noting that 


t ^ B 2 u |'-v T j»-mB 2 (u T .,-V|i) 

l v z'h. = 


~ 2 -^ 2 


-m B 


and along the leading edge for an infinitesimal strip of width At]’', 

U; , = V. = 0 

t i ’ 
u o 

V* 1 ’ ■ T 


V* 1 ’ - - u o 


2121 
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results in 


u 0 Vl-B 2 m 2 
= K - 



r’ 


( 7 ) 


At subsonic speeds, the upper limit = 00 , and the integral is 

divergent. In actual cases, however, this divergence does not present 
any diff iculty because the construction of the finite leading edge by 
means of the superposition of two semi-infinite leading edges of 
opposite sign leads to the result i that the infinite upper limit 
cancels. At supersonic speeds, Sg is the position of the last 
source with Mach after cone including the point x',y',z'j that is, 

£g = x'-B Vy ,2 +z' 2 

Integration of equation (7) yields the following expressions for 
the line source originating at (0,0): 

At subsonic speeds, 


v T (x,y,z) = 


u 0 Vl +in2 Bi 2 


2jtm 


— 1 X * 

sinh x i- 
5 A 


(8a) 


where 

v d = VI B2 l(y ,2+z ' 2 ) 


and 



In equation (8a), the term that arises from substituting the upper 
limit of integration has been neglected because, as indicated pre- 
viously, it vanishes for leading edges of finite lengths. 

At supersonic speeds, for Bm < 1 (subsonic leading edge), 

w T (x,y,z) = — Vl-B^m 2 cosh -1 ri (8b) 
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and for Bm > 1 (supersonic leading edge), equation (8b) becomes 

w L (x,y,z) = - ~ V® 2 ™ 2 - 1 cos" 1 FT (8c) 

a 

In equations (8) and in all subsequent expressions, the positive 
value in a radical term must be preserved when extracting the root. 

For example, if 

y < 0 

then 1 

= V( -y ) 2 = -y = |y| 

For a leading edge normal to the flight direction, m = » and 
equations (8) yield the following expressions: 

At subsonic speeds, 

B i u n , B^y 

w L (x,y,z) = — sinh -1 - j - (8d) 

At supersonic speeds , 

v L (x,y,z) -2 cos' 1 ~ By (8e) 

Vx 2 -B 2 s 2 

It is shown by equations (8) that the downvash fields w^ con- 
tributed by the oblique and normal leading edges are conical, or w 
is constant along radial lines emanating from the origin. - 

Contribution of leading-edge discontinuity at origin . - From fig- - 
ure 2, it can be seen that the leading edge cut-off at the origin 
results in 

V--T (9a) 

along the semi-infinite side edge. When equation (9a) is substituted 
into equation (3) (with u, equal to zero along the side edge), there 
results ’ ' ■ ■ * 
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fz 

w LD (x,y,z) = - (9b) 

where t 2 is the position of the last source that can influence the 
point x,y,z; that is, at subsonic speeds, 4 2 i- s & t infinity; and at 
supersonic speeds , 

1 2 = x-B Vy 2 +z 2 


Integration of equation (9b) yields the following: 

i 

At subsonic speeds, 

v LD ( X< y,2) = - sinh' 1 £ (10a) 

where the term arising from the infinite upper limit of integration has 
been neglected because it vanishes for leading edges of finite length 
and 

e d = V|b 2 | (y 2 +z 2 ) 

At supersonic speeds, 


u. 


v T _(x,y,z) cosh -1 

LD ran & d 


(10b) 


It is shown by equations (10) that the downwash contribution of the 
leading-edge discontinuity at the origin results entirely from the 
obliquity of the leading edge.. Thus, if the leading edge is normal to 
the flight direction, this contribution vanishes. 


The Semi-Infinite Streamwise Side Edge 

Along the streamwise side edge, the quantity u. vanishes; thus 
the source-distribution function along this edge is proportional only ' 
to the quantity v^. 

' The details of evaluating the contribution of the streamwise side 
edge to the downwash field are given in the appendix. The results for 
the semi-infinite left streamwise side edge are as follows : , 



12 


NACA TN 2344 


At subsonic speeds, 


w s (x,y,z) 


hS (/ f ' 1 2x_\ 

Zn \x + Vx 2 +I d 2 la / 


(Ha) 


At supersonic speeds. 


w s (x,y,z) 


B 2 u 0 y f l 
51 yx + Vx 2 -ld 2 



(lib) 


Equations (11) show that the dovnvash contribution of the strearawise 
side edge is independent of the obliquity of the leading edge. 


FINITE WINGS WITH UNIFORM LOADING 

The dovnvash field due to finite wings with uniform loading can 
be obtained by superimposing the fields due to a number of semi-infinite 
wings of the type considered in the previous section (fig. l) . This 
superposition is equivalent to superimposing the fields due to the 
source lines expressed by equations (8), (10), and (11). 

A plan form with curved edges requires an infinite number of 
source lines. If the edges are composed of straight-line segments, 
however, a finite number of source lines can be used to represent the 
plan form. 


Swept Wings with Streamwise Tips 

Let the downwash field due to a semi-infinite line source 
beginning at x n ,y n be denoted by w ■ s let the subscripts 

0 0 x o ,y o 

L, S, T, LD, and TD refer to the leading, side, and trailing edges 
and the leading- and trailing-edge discontinuities, respectively, all for 
the right half -wing; and let w refer to the effect of the source 
lines originating on the left half -wing. Then, the dovnvash field for 
a uniformly loaded wing, which is symmetrical with respect to the 
x-axis and which has uniformly swept leading edges and streamwise tips, 
is given by the following sum: 
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In equation (12a) and in the subsequent expressions, the upper and 
lower signs preceding a term refer to the right and left half -wings, 
respectively. 

The downwash fields due to the semi -infinite line sources indi- 
cated in equation (12a) can be obtained by simple manipulation of equa- 
tions (8), (10), and (ll) . The following transformations are made to 
obtain the effect of a semi-infinite line source originating at Xq^q: 


in equations (8), (10) , and (11) = 


in equations (8), (10), and (ll) = -(y-yo) 


Equations (8) and (10) can also be applied to a semi-infinite line 
source along the trailing edge by replacing m in these equations 
with the cotangent of the sweep angle of ' the trailing edge. 

At subsonic speeds, every point in the field is affected by all 
the terms in equation (12a); whereas at supersonic speeds, the point 
is affected only by those terms that refer to edges which lie within 
the Mach forecone from the point under consideration. 


Rectangular Wings 

For the rectangular wing, the edge discontinuities LD and TD 
disappear and equation (12a) becomes 
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WINGS WITH ARBITRARY LOADING 

If the wing loading is an arbitrary function of x and y, the 
integration in equation (3) for the downwash field is, in general., 
required over the entire wing surface and in the wake. An alternative 
•procedure in this case is the use of "correspondence formulas" as 
indicated in reference 8 by means of which the downwash field due to 
a variable loading may be expressed in terms of the downwash field due 
to a uniform loading with the addition, of corrections for the edges of 
the plan form. 

In reference 8, it is shown that if u on the wing is expressed 
in a series as a function of x and y in terms of the uniform 
prescribed velocity uq and any term in the series is differentiated 
with respect to x and y, s and t times, .respectively, such that 

(u) s t = bu o (° n w *- n s) ( 13a ) 

x y 

where b is a constant; then subject to edge corrections, there 
results 

( w n,k) s t ( X >Y> Z ) = bw o( x ' y ' z ) ( 13b ) 

x y 

where w , and w^ refer to the downwash fields due to the term in 

iij k 0 

the series and to a uniform loading, respectively. The edges of finite 
plan forms may alter the given relation (13a) so that edge corrections 
may be required for the relation expressed by equation (13b). At 
supersonic speeds, the relation given by equation (13b) thus applies 
to finite wings at all points in the flow field outside of the Mach 
aftercone from the edges that alter the given relation (13a) oh the 
wing. Correspondence formulas for rectangular wings at subsonic and 
supersonic speeds are given in reference 8, table I. 
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ILLUSTRATIVE EXAMPLE FOR RECTANGULAR WING 
AT SUPERSONIC SPEEDS 

As an illustration of the method described herein, the downwash 
field is obtained at supersonic speeds for a rectangular wing with 
uniform loading and with linear chordwise variation in loading. The 
lettered regions in the subsequent discussion refer to figure 3. 


Uniform Loading 

For the rectangular wing with uniform loading at supersonic speed, 
equations (8e), (lib), and (12c) are applicable. 

Region A . - In the region within the leading-edge Mach cones and 
outside of the side-edge and trailing-edge Mach cones, 

i* 0 )iS x >?> z ) = ( W 0,0 + *0,0>L 

= - /cos" 1 ~~ ■ - + COS -1 By ^ = -Bu 0 (14) 

51 \ BV 

Region B . - In the region within the trailing-edge Mach cones and 
outside of the side -edge Mach cones, 

(v 0 ) B (x,y,z) = (w 0 )a " (w c ,0 + *c,oh 

By ) =0 
x- c) 2 -B 2 z 2 / 

(15) 

Region C . - Region C, which refers to points outside the trailing- 
edge Mach cones, is divided into two subregions depending on whether 
the point lies within one or both side-edge Mach cones. 

Region C-j_. - If the point lies within one side-edge Mach cone. 


= -Bu,. 


Bu 0 / 

+ cos 


-By 


* V 


2 _2 2 
B z 


+ cos 


-1 
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(v 0 ) C i(x,y,z) 


(w 0 ) A - (v 0 , h ) L 


( w O,h)s 


Mh + !?o 

'"O' A ^ Jt 


-1 ~Bya 

COS a/'§ 'T ' 2 " By « 




Bu, 


0 


_-l By a 
cos — — +i3y. 


Vx 2 -B 2 z^ \x+V x2 -^ a 2 


(16) 


where 


S a = B^y^+z 2 


and 


y a = y-h 

Region Cg. - For points within both side-edge Mach cones, 


(v 0 ) C2 (x,y,z) • (v^ - (w 0> ; h ) L -.('' 0j . h > s 


Bur 


“ (Vc! + n 


_i Byb 

cos ~7rT = T9 +By B 
Vx 2 - B d z^ 


l x+Vx 2 -U 2 £y > 2 


where 


x 

b s b* 

(17) 


e b = b Vyb 2+z2 


and 


y b = y+h 

Region D. - Region D refers to points that are always within Mach 
cones from, the leading, side, and trailing edges. This region is sub- 
divided into five regions, as shown in figure 3. 


Region Dj_. - For points in 

, (v 0 ) Di (x,y,z) = (v Q ) Ci - (w c>Q + w c>0 ) T = (w Q ) Ci + Bu q (18) 
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Region D2. - For points in Dg> 


('' 0 ) D s (x ’ y ’ 2) = ("o'c, + + tB c,h>T + (v c,h>S 


/x bu 0 -1 fi y a 

(w V + —— cos — 

0 C 1 ” A/(x-0 2 -B 2 2 2 


x-c+A/(x-c) 2 -£ a 2 S a 


Region D3. - For points in D3, 
( w 0^D 3 ( Xy ^ z ) - - <* 0 ,-h>L ’ ^0,-^S 


( ) + £2 cos- 1 

0Dl ” V^V 


x+ Vx - 


2 S 2 


Region D 4 . - For points in D 4 , 


(w 0 )D 4 (x,y,z) = (w 0 ) D2 - (v 0 ^_ h ) L - (w 0 ^. h ) s = (v 0 ) D2 + (wq)^ - (vq)^ 


Region D5. - For points in D5, 


(w 0 ) D (x,y,z) = (w 0 ) D + (w c ^. h ) T + (w c ^ h ) s 


/ \ Bu 0 -1 
( w o)d 4 - — cos 


/(x- c) 2 -B 2 z 2 


x-c+A/(x-c) 2 -e b 2 ^ 
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Linear Chordwise Loading 
For linear chordwise loading, let 

= au Q x (on wing) (23a) 


where a and Uq are constants. The downwash field for this type of 
loading can be obtained by using the correspondence formulas of refer- 
ence 8 in conjunction with the downwash field obtained in the preceding 
section for the uniformly loaded rectangular wing at supersonic speed. 
If w^ and Wq refer to the downwash fields due to linear chordwise 

loading and to uniform loading, respectively, then from reference 8, 
table I, for all points outside the trailing-edge Mach cones, 

v x (x,y,z) => a f v 0 (£,y,z) d£ (23b) 

Jbz 


and for points within the trailing-edge Mach cones. 


w x (x,y,z) 



woU>y>z) d£+ c(Aw Q ) T (x,y,z) 


(23c) 


where (Awq)^ refers to the effect of the wing cut-off at the trailing 
edge for uniform loading. 

Region A . - For the region within the leading-edge Mach cone and 
outside the side-edge and trailing-edge Mach cones, equations (23b) 
and (14) are applicable and there results 

(wi) A (x>y>z) = -aBu 0 (x-Bz). (24) 

Region B . - In the region within the trailing-edge Mach cones and 
outside the side-edge Mach cones, equations (23c), (14), and (15) yield 

J nc+Bz 

d£ + acBu 0 = 0 (25) 

Bz 

Region C . - In the region within one or both side-edge Mach cones 
and outside the trailing-edge Mach cones, equation (23b) applies. 
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Region C^. - For points in region C-^, equation (23b) may be 
written 


(vi) Cl (x,y,z) 


~ a <f w o (e ^ z) 



<u 



(26a) 


In equation (26a), w 0 (£,y,z) is represented as (w Q ) Cl (equation (16)) 
throughout the entire range of integration because (wq) ^ evaluated 
in region A is equal to (wq)^ (equation (14)), inasmuch as the 
imaginary part is discarded. 


Equation (26a), requires the evaluation of the following integrals: 

(26b) 


«Bz 


-1 ®y a . t 
cos — ■ d£ 


Vi 2 - B 2 z 2 


j2„2 


where By a and Bz are constants, and 


12 


*l a 


e £ 


di 


(26c) 


where i a is a constant. 


The solution for 1^ obtained by integration by parts and simpli- 
fication is 
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1 B y a lx 1 By a x 

I(x,y ,z)=xcos' _____ -By cosh - j~ - Bz cos" — = 

1 Vx 2 -B 2 z 2 U Vx 2 -B 2 z 2 


The solution for I 2 is 


I 2 ( x >y a ^ z ) - 


cosh' 


•1 x 


\ l 


V x2 - 


7 

a _J 


(26d) 


(26e) 


Utilizing the solutions for Ij and I 2 given by equations (26d) 
and (26e), equation (26a) yields for points in region C within one side- 
edge Mach cone, 

a ® u o r n 

(Vi) Cl (x,y^) « - If — Lli(x,y a ,z) + By a I 2 (x,y a ,z) I (26f) 

Region C 2 . - For points in region C within both side-edge Mach 
cones, equation (23b) may be written 

>x 


( w lW x >y> z ) = a / w (e,y,z) d £ 
* J Bz 

7 

On? 


= ( w i) c x + — 


-1 By b 

cos — - d£ + 


V^ 2 -b 


2,2 


•' C fcvm? ■ 5 s ) 


(27a) 


where wq(£ , y, z) has been obtained from equations (16) and (17). 
Utilizing the solutions for 1^. and I 2 given in equations (26d) 
and (26e), in which y and £ a are replaced by y b and £ b , 
respectively, changes equation (27a) to yield 

aBu 0 r -i 

( w 1 )c 2 (x^y>z) = (w x ) Ci + ~ y ~ [_Il(x,y b ,z) + By^ I 2 (x,y b ,z)J 


(27b) 
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Region D . - In the region within the trailing-edge Mach cone and 
within one or both side-edge Mach cones, equation (23c) applies. 

Region D^. - For points in region D^, the integral in equa- 
tion (23c) becomes 

J ^x pX 

W Q ( d£ = (w 1 ) c + aBu Q I di- = (w 1 ) + aBu (x-c-Bz) 

Bz 1 Jc+Bz 1 


( 2£a) 

where Wq(4,y,z) has been obtained from equations (16) and (18). 

The quantity (Awq)^ in equation (23c) is evaluated for this case 
as the 'difference between ( w q)d-^ and (wq) that is, from equa- 
tion (18), (AWq)^, = BUq. When these results are combined, equa- 
tion (23c) yields 


(wi) Dl (x,y,z) = (w^^ + aBu Q (x-Bz) 


Region D 2 . - For points in region D 2 , the integral in equa- 
tion (23c) becomes 


JV 

«JBz 


aBu 0 

y,z)dt= (v x ) Ci + — 


f ■ 

*-*c+Bz 


By a 


^S-c) 2 - 


B 2 z 2 


u + 


(>28b) 


By fl 


L+i \ ^- c+ V ^“ c ) 2 '^ a 2 



de 


(29a) 


where w Q ( £,y,z) has been obtained from equations (16) and (19). The 

integrals in equation (29a) may be evaluated by utilizing the solutions 
for Ij_ and I 2 as expressed by equations (26d) and (26e). In these 
solutions x is replaced by x-c. The quantity (Aw 0 ) t in equa- 
tion (23c) is evaluated in this case as the difference between (w ) 

0 

and (w Q ) c . When these results are combined, equation (23c) yields 
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aBu, 


(v^d/x^z) = (w x ) Ci + — 


Il(x-c,y a ,z) + By a I 2 (x-c,y a ,z) + 


cos 


-1 


By a 


+By. 


x-c 


V(x- c) 2 -B 2 z 2 a \x-c+V(x-c) 2 -e a ^ ^71 


(29b) 


Region D 3 . - For points in region D 3 , equation (23c) yields 


(vi) D3 (x,y,z) = (w 1 ) Di + 


By b 


Y 

/ COE 

J Bz 

X> (^¥V 


7 e 2 -B 2 z 2 

j; 


ae + 

d£ 


(30a) 


where w Q (E,y,z) has been obtained from equations (18) and (20) and the 
term (Aw 0 ) t is included in (w 1 ) D . The integrals in equation (30a) 

may be evaluated by utilizing the solutions for I-, and Ig as 
expressed by equations (26d) and (26e). In these solutions, y a is 
replaced by y b ; thus 

aBUp r — * 

( w i)d 3 (^^ z ) = ( w i) Dl + ~T~ |_ I l (x ^ y b ,z) + By b i 2 (x ^ y b^ z )J (30b) 
Region D 4 . - For points in region D 4 , equation (23c) yields 

<"i) D4 U,y,2) * < u i>d 2 + < v i>d 3 - (v 1 , d 1 

where w Q (i,y,z) has been obtained from equations (19) and (21). 


( 31 ) 
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Region Dg. - For points in region Dg, equation (23c) yields 


aBu r 


(vi) D g(x,y, Z ) - (w^ ~ I 008 

*c+Bz 


r 


-i B yb 


ia -c) 2 -b 2 z 2 


ai + 


By* 


r r • i 

J c+ u v- c+/ V ( M 2 -e b 2 ^ 


d£ + 


cos 


-1 


Byb 


Y(x-c) 2 -B 2 z 2 


+ Byb 


x-c+ 


V(x- C ) 2 -e b 2 ^ 


x-c 

2 j 


(32a) 

where WQ(£,y,z) has been obtained from equation (22). The integrals 
in equation (32a) may be evaluated by utilizing the solutions for 1^ and. 
1 2 as expressed by equations (26d) and (26e). In these solutions, 
x is replaced by x-c, and y a and £ a are replaced by y b and £ b , 
respectively. Thus 

aBug 

( w l) D g(x,y,z) = (w-l)^ - -j— ^ I 1 (x-c,y b ,z) + By b I 2 (x-c,y b ,z) + 


COS 


-1 5& 


~*' D + gy | ± X-c | I V 

V(x- c) 2 -B 2 Z 2 \X-C+^/(x-c ) 2 -^ b 2 y/jj 


(32b) 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 22, 1951. 
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APPENDIX 

DERIVATION OF DOWNWASH FIELD DUE TO SEMI-INFINITE 
LEFT STREAMWISE SIDE EDGE 


The subsequent derivation for the downwash field refers to the 
left streamwise side edge of the uniformly loaded semi-infinite oblique 
wing, as shown in figure 1. 

The downwash field contributed by the side edge may be obtained by 
means of equation (3), in which the integral is evaluated along the side 
edge. As noted in the text, the quantity U£ is zero along the stream- 
wise side edge (of width At))j therefore. 


w g (x,y,z) 



• (Al) 


Across 
increase is 
side edge. 


the side edge there is an abrupt increase in 
assumed to occur over an infinitesimal width 
Then according to equation (4a), 

■»x 


v(x,0,0) 


-i; 


VLr 


U, 


&t= -j- 


u. This ' 
e across the 


(A2) 


At 


y = £ , however, u^ 


due to the side edge vanishes and therefore 


v(x,c,0) = 0 


(A3) 


The distribution- of v in the vicinity of the side edge (y = 0) 
is shown in the following sketch: 



y 3** 
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It is assumed that the abrupt changes in v indicated in the 
sketch occur over the infinitesimal distances An, so that along the 
line (I >0,0), 


= 


u o e 


(A4) 


and along the line (£,e,0). 


„ u 0 e 
V> r 


(AS) 


Substituting these values for v^A*! into equation (Al) results in 


Ku^ 

w s (x,y,z) = — 


1 


x-B Vy^+z^ 


i 


V(x4) 2 - B 2 (y 2 +z 2 ) 


f 

do 


c, x-B V(y-e) 2 +z^ 


iai 

V(x4) 2 -B 2 (y-e) 2 +z 2 


(A6) 


where the two upper limits refer to subsonic and supersonic speeds, 
respectively. Integration of equation (A6) yields 


w g (x,y,z) 


Ku 

T 


^,x-B 


— \^J X log I X "^ + AA x "^ ) 2 -B 2 (y 2 +z 2 ) | + 

V(x4) 2 - B 2 (y 2 +z 2 } 

J 0 

{x log | x-£ + a/( x-0 2 -B 2 [( y-€ ) 2 +z^] | 
j*>,x-B V( y-e) 2 +z 2 

V( x-0 2 -B 2 [( y-e ) 2 +z3 [ 


(A7) 
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When the limits are substituted into equation (A7) and e is 
made to approach zero, the following results are obtained: 

At subsonic speeds. 


•w s (x,y,z) 


B i 2u py 

2n 



where 


* V|b 2 | (y 2 +z 2 ) 


(A8) 


At supersonic speeds, 


w s (x,y,z) 



(A9) 
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PiKure 3. - Plan view of regions considered in examples for rectangular wing at supersonic speeds. Regions bounded by portions of three-dimensional 

& Ma aft*. (<nnpa. 



